Improved metagenome assemblies through selective enrichment of bacterial genomic DNA from eukaryotic host genomic DNA using ATAC-seq.

Genomics can be used to study the complex relationships between hosts and their microbiota. Many bacteria cannot be cultured in the laboratory, making it difficult to obtain adequate amounts of bacterial DNA and to limit host DNA contamination for the construction of metagenome-assembled genomes (MAGs). For example, Wolbachia is a genus of exclusively obligate intracellular bacteria that live in a wide range of arthropods and some nematodes. While Wolbachia endosymbionts are frequently described as facultative reproductive parasites in arthropods, the bacteria are obligate mutualistic endosymbionts of filarial worms. Here, we achieve 50-fold enrichment of bacterial sequences using ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) with Brugia malayi nematodes, containing Wolbachia (wBm). ATAC-seq uses the Tn5 transposase to cut and attach Illumina sequencing adapters to accessible DNA lacking histones, typically thought to be open chromatin. Bacterial and mitochondrial DNA in the lysates are also cut preferentially since they lack histones, leading to the enrichment of these sequences. The benefits of this include minimal tissue input (<1 mg of tissue), a quick protocol (<4 h), low sequencing costs, less bias, correct assembly of lateral gene transfers and no prior sequence knowledge required. We assembled the wBm genome with as few as 1 million Illumina short paired-end reads with >97% coverage of the published genome, compared to only 12% coverage with the standard gDNA libraries. We found significant bacterial sequence enrichment that facilitated genome assembly in previously published ATAC-seq data sets from human cells infected with Mycobacterium tuberculosis and C. elegans contaminated with their food source, the OP50 strain of E. coli. These results demonstrate the feasibility and benefits of using ATAC-seq to easily obtain bacterial genomes to aid in symbiosis, infectious disease, and microbiome research.

Common analysis of direct RNA sequencinG CUrrently leads to misidentification of m5C at GCU motifs.

RNA modifications, such as methylation, can be detected with Oxford Nanopore Technologies direct RNA sequencing. One commonly used tool for detecting 5-methylcytosine (m5C) modifications is Tombo, which uses an "Alternative Model" to detect putative modifications from a single sample. We examined direct RNA sequencing data from diverse taxa including viruses, bacteria, fungi, and animals. The algorithm consistently identified a m5C at the central position of a GCU motif. However, it also identified a m5C in the same motif in fully unmodified in vitro transcribed RNA, suggesting that this is a frequent false prediction. In the absence of further validation, several published predictions of m5C in a GCU context should be reconsidered, including those from human coronavirus and human cerebral organoid samples.

A novel family of sugar-specific phosphodiesterases that remove zwitterionic modifications of GlcNAc.

The zwitterions phosphorylcholine (PC) and phosphoethanolamine (PE) are often found esterified to certain sugars in polysaccharides and glycoconjugates in a wide range of biological species. One such modification involves PC attachment to the 6-carbon of N-acetylglucosamine (GlcNAc-6-PC) in N-glycans and glycosphingolipids (GSLs) of parasitic nematodes, a modification that helps the parasite evade host immunity. Knowledge of enzymes involved in the synthesis and degradation of PC and PE modifications is limited. More detailed studies on such enzymes would contribute to a better understanding of the function of PC modifications and have potential application in the structural analysis of zwitterion-modified glycans. In this study, we used functional metagenomic screening to identify phosphodiesterases encoded in a human fecal DNA fosmid library that remove PC from GlcNAc-6-PC. A novel bacterial phosphodiesterase was identified and biochemically characterized. This enzyme (termed GlcNAc-PDase) shows remarkable substrate preference for GlcNAc-6-PC and GlcNAc-6-PE, with little or no activity on other zwitterion-modified hexoses. The identified GlcNAc-PDase protein sequence is a member of the large endonuclease/exonuclease/phosphatase superfamily where it defines a distinct subfamily of related sequences of previously unknown function, mostly from Clostridium bacteria species. Finally, we demonstrate use of GlcNAc-PDase to confirm the presence of GlcNAc-6-PC in N-glycans and GSLs of the parasitic nematode Brugia malayi in a glycoanalytical workflow.

Defining the filarial N-glycoproteome by glycosite mapping in the human parasitic nematode Brugia malayi.

N-linked glycosylation is a critical post translational modification of eukaryotic proteins. N-linked glycans are present on surface and secreted filarial proteins that play a role in host parasite interactions. Examples of glycosylated Brugia malayi proteins have been previously identified but there has not been a systematic study of the N-linked glycoproteome of this or any other filarial parasite. In this study, we applied an enhanced N-glyco FASP protocol using an engineered carbohydrate-binding protein, Fbs1, to enrich N-glycosylated peptides for analysis by LC-MS/MS. We then mapped the N-glycosites on proteins from three host stages of the parasite: adult female, adult male and microfilariae. Fbs1 enrichment of N-glycosylated peptides enhanced the identification of N-glycosites. Our data identified 582 N-linked glycoproteins with 1273 N-glycosites. Gene ontology and cell localization prediction of the identified N-glycoproteins indicated that they were mostly membrane and extracellular proteins. Comparing results from adult female worms, adult male worms, and microfilariae, we find variability in N-glycosylation at the protein level as well as at the individual N-glycosite level. These variations are highlighted in cuticle N-glycoproteins and adult worm restricted N-glycoproteins as examples of proteins at the host parasite interface that are well positioned as potential therapeutic targets or biomarkers.

Unraveling cross-reactivity of anti-glycan IgG responses in filarial nematode infections.

Parasitic nematodes responsible for filarial diseases cause chronic disablement in humans worldwide. Elimination programs have substantially reduced the rate of infection in certain areas, but limitations of current diagnostics for population surveillance have been pointed out and improved assays are needed to reach the elimination targets. While serological tests detecting antibodies to parasite antigens are convenient tools, those currently available are compromised by the occurrence of antibodies cross-reactive between nematodes, as well as by the presence of residual antibodies in sera years after treatment and clearance of the infection. We recently characterized the N-linked and glycosphingolipid derived glycans of the parasitic nematode Brugia malayi and revealed the presence of various antigenic structures that triggered immunoglobulin G (IgG) responses in infected individuals. To address the specificity of IgG binding to these glycan antigens, we screened microarrays containing Brugia malayi glycans with plasma from uninfected individuals and from individuals infected with Loa loa, Onchocerca volvulus, Mansonella perstans and Wuchereria bancrofti, four closely related filarial nematodes. IgG to a restricted subset of cross-reactive glycans was observed in infection plasmas from all four species. In plasma from Onchocerca volvulus and Mansonella perstans infected individuals, IgG binding to many more glycans was additionally detected, resulting in total IgG responses similar to the ones of Brugia malayi infected individuals. For these infection groups, Brugia malayi, Onchocerca volvulus and Mansonella perstans, we further studied the different IgG subclasses to Brugia malayi glycans. In all three infections, IgG1 and IgG2 appeared to be the major subclasses involved in response to glycan antigens. Interestingly, in Brugia malayi infected individuals, we observed a marked reduction in particular in IgG2 to parasite glycans post-treatment with anthelminthic, suggesting a promising potential for diagnostic applications. Thus, we compared the IgG response to a broad repertoire of Brugia malayi glycans in individuals infected with various filarial nematodes. We identified broadly cross-reactive and more specific glycan targets, extending the currently scarce knowledge of filarial nematode glycosylation and host anti-glycan antibody response. We believe that our initial findings could be further exploited to develop disease-specific diagnostics as part of an integrated approach for filarial disease control.

Alteration of rhesus macaque serum N-glycome during infection with the human parasitic filarial nematode Brugia malayi.

Serum N-glycan profiling studies during the past decades have shown robust associations between N-glycan changes and various biological conditions, including infections, in humans. Similar studies are scarcer for other mammals, despite the tremendous potential of serum N-glycans as biomarkers for infectious diseases in animal models of human disease and in the veterinary context. To expand the knowledge of serum N-glycan profiles in important mammalian model systems, in this study, we combined MALDI-TOF-MS analysis and HILIC-UPLC profiling of released N-glycans together with glycosidase treatments to characterize the glycan structures present in rhesus macaque serum. We used this baseline to monitor changes in serum N-glycans during infection with Brugia malayi, a parasitic nematode of humans responsible for lymphatic filariasis, in a longitudinal cohort of infected rhesus macaques. Alterations of the HILIC-UPLC profile, notably of abundant structures, became evident as early as 5 weeks post-infection. Given its prominent role in the immune response, contribution of immunoglobulin G to serum N-glycans was investigated. Finally, comparison with similar N-glycan profiling performed during infection with the dog heartworm Dirofilaria immitis suggests that many changes observed in rhesus macaque serum N-glycans are specific for lymphatic filariasis.

Mass Spectrometric and Glycan Microarray-Based Characterization of the Filarial Nematode Brugia malayi Glycome Reveals Anionic and Zwitterionic Glycan Antigens.

Millions of people worldwide are infected with filarial nematodes, responsible for lymphatic filariasis (LF) and other diseases causing chronic disablement. Elimination programs have resulted in a substantial reduction of the rate of infection in certain areas creating a need for improved diagnostic tools to establish robust population surveillance and avoid LF resurgence. Glycans from parasitic helminths are emerging as potential antigens for use in diagnostic assays. However, despite its crucial role in host-parasite interactions, filarial glycosylation is still largely, structurally, and functionally uncharacterized. Therefore, we investigated the glycan repertoire of the filarial nematode Brugia malayi. Glycosphingolipid and N-linked glycans were extracted from several life-stages using enzymatic release and characterized using a combination of MALDI-TOF-MS and glycan sequencing techniques. Next, glycans were purified by HPLC and printed onto microarrays to assess the host anti-glycan antibody response. Comprehensive glycomic analysis of B. malayi revealed the presence of several putative antigenic motifs such as phosphorylcholine and terminal glucuronic acid. Glycan microarray screening showed a recognition of most B. malayi glycans by immunoglobulins from rhesus macaques at different time points after infection, which permitted the characterization of the dynamics of anti-glycan immunoglobulin G and M during the establishment of brugian filariasis. A significant level of IgG binding to the parasite glycans was also detected in infected human plasma, while IgG binding to glycans decreased after anthelmintic treatment. Altogether, our work identifies B. malayi glycan antigens and reveals antibody responses from the host that could be exploited as potential markers for LF.

X-treme loss of sequence diversity linked to neo-X chromosomes in filarial nematodes.

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

Dual RNAseq analyses at soma and germline levels reveal evolutionary innovations in the elephantiasis-agent Brugia malayi, and adaptation of its Wolbachia endosymbionts.

Brugia malayi is a human filarial nematode responsible for elephantiasis, a debilitating condition that is part of a broader spectrum of diseases called filariasis, including lymphatic filariasis and river blindness. Almost all filarial nematode species infecting humans live in mutualism with Wolbachia endosymbionts, present in somatic hypodermal tissues but also in the female germline which ensures their vertical transmission to the nematode progeny. These α-proteobacteria potentially provision their host with essential metabolites and protect the parasite against the vertebrate immune response. In the absence of Wolbachia wBm, B. malayi females become sterile, and the filarial nematode lifespan is greatly reduced. In order to better comprehend this symbiosis, we investigated the adaptation of wBm to the host nematode soma and germline, and we characterized these cellular environments to highlight their specificities. Dual RNAseq experiments were performed at the tissue-specific and ovarian developmental stage levels, reaching the resolution of the germline mitotic proliferation and meiotic differentiation stages. We found that most wBm genes, including putative effectors, are not differentially regulated between infected tissues. However, two wBm genes involved in stress responses are upregulated in the hypodermal chords compared to the germline, indicating that this somatic tissue represents a harsh environment to which wBm have adapted. A comparison of the B. malayi and C. elegans germline transcriptomes reveals a poor conservation of genes involved in the production of oocytes, with the filarial germline proliferative zone relying on a majority of genes absent from C. elegans. The first orthology map of the B. malayi genome presented here, together with tissue-specific expression enrichment analyses, indicate that the early steps of oogenesis are a developmental process involving genes specific to filarial nematodes, that likely result from evolutionary innovations supporting the filarial parasitic lifestyle.

The Wolbachia Symbiont: Here, There and Everywhere.

Wolbachia symbionts, first observed in the 1920s, are now known to be present in about 30-70% of tested arthropod species, in about half of tested filarial nematodes (including the majority of human filarial nematodes), and some plant-parasitic nematodes. In arthropods, they are generally viewed as parasites while in nematodes they appear to be mutualists although this demarcation is not absolute. Their presence in arthropods generally leads to reproductive anomalies, while in nematodes, they are generally required for worm development and reproduction. In mosquitos, Wolbachia inhibit RNA viral infections, leading to populational reductions in human RNA virus pathogens, whereas in filarial nematodes, their requirement for worm fertility and survival has been channeled into their use as drug targets for filariasis control. While much more research on these ubiquitous symbionts is needed, they are viewed as playing significant roles in biological processes, ranging from arthropod speciation to human health.

Diminutive, degraded but dissimilar: Wolbachia genomes from filarial nematodes do not conform to a single paradigm.

Wolbachia are alpha-proteobacteria symbionts infecting a large range of arthropod species and two different families of nematodes. Interestingly, these endosymbionts are able to induce diverse phenotypes in their hosts: they are reproductive parasites within many arthropods, nutritional mutualists within some insects and obligate mutualists within their filarial nematode hosts. Defining Wolbachia 'species' is controversial and so they are commonly classified into 17 different phylogenetic lineages, termed supergroups, named A-F, H-Q and S. However, available genomic data remain limited and not representative of the full Wolbachia diversity; indeed, of the 24 complete genomes and 55 draft genomes of Wolbachia available to date, 84 % belong to supergroups A and B, exclusively composed of Wolbachia from arthropods. For the current study, we took advantage of a recently developed DNA-enrichment method to produce four complete genomes and two draft genomes of Wolbachia from filarial nematodes. Two complete genomes, wCtub and wDcau, are the smallest Wolbachia genomes sequenced to date (863 988 bp and 863 427 bp, respectively), as well as the first genomes representing supergroup J. These genomes confirm the validity of this supergroup, a controversial clade due to weaknesses of the multilocus sequence typing approach. We also produced the first draft Wolbachia genome from a supergroup F filarial nematode representative (wMhie), two genomes from supergroup D (wLsig and wLbra) and the complete genome of wDimm from supergroup C. Our new data confirm the paradigm of smaller Wolbachia genomes from filarial nematodes containing low levels of transposable elements and the absence of intact bacteriophage sequences, unlike many Wolbachia from arthropods, where both are more abundant. However, we observe differences among the Wolbachia genomes from filarial nematodes: no global co-evolutionary pattern, strong synteny between supergroup C and supergroup J Wolbachia, and more transposable elements observed in supergroup D Wolbachia compared to the other supergroups. Metabolic pathway analysis indicates several highly conserved pathways (haem and nucleotide biosynthesis, for example) as opposed to more variable pathways, such as vitamin B biosynthesis, which might be specific to certain host-symbiont associations. Overall, there appears to be no single Wolbachia-filarial nematode pattern of co-evolution or symbiotic relationship.

Population genomic evidence that human and animal infections in Africa come from the same populations of Dracunculus medinensis.

BACKGROUND: Guinea worm-Dracunculus medinensis-was historically one of the major parasites of humans and has been known since antiquity. Now, Guinea worm is on the brink of eradication, as efforts to interrupt transmission have reduced the annual burden of disease from millions of infections per year in the 1980s to only 54 human cases reported globally in 2019. Despite the enormous success of eradication efforts to date, one complication has arisen. Over the last few years, hundreds of dogs have been found infected with this previously apparently anthroponotic parasite, almost all in Chad. Moreover, the relative numbers of infections in humans and dogs suggests that dogs are currently the principal reservoir on infection and key to maintaining transmission in that country. PRINCIPAL FINDINGS: In an effort to shed light on this peculiar epidemiology of Guinea worm in Chad, we have sequenced and compared the genomes of worms from dog, human and other animal infections. Confirming previous work with other molecular markers, we show that all of these worms are D. medinensis, and that the same population of worms are causing both infections, can confirm the suspected transmission between host species and detect signs of a population bottleneck due to the eradication efforts. The diversity of worms in Chad appears to exclude the possibility that there were no, or very few, worms present in the country during a 10-year absence of reported cases. CONCLUSIONS: This work reinforces the importance of adequate surveillance of both human and dog populations in the Guinea worm eradication campaign and suggests that control programs aiming to interrupt disease transmission should stay aware of the possible emergence of unusual epidemiology as pathogens approach elimination.

A Meta-Analysis of Wolbachia Transcriptomics Reveals a Stage-Specific Wolbachia Transcriptional Response Shared Across Different Hosts.

Wolbachia is a genus containing obligate, intracellular endosymbionts with arthropod and nematode hosts. Numerous studies have identified differentially expressed transcripts in Wolbachia endosymbionts that potentially inform the biological interplay between these endosymbionts and their hosts, albeit with discordant results. Here, we re-analyze previously published Wolbachia RNA-Seq transcriptomics data sets using a single workflow consisting of the most up-to-date algorithms and techniques, with the aim of identifying trends or patterns in the pan-Wolbachia transcriptional response. We find that data from one of the early studies in filarial nematodes did not allow for robust conclusions about Wolbachia differential expression with these methods, suggesting the original interpretations should be reconsidered. Across datasets analyzed with this unified workflow, there is a general lack of global gene regulation with the exception of a weak transcriptional response resulting in the upregulation of ribosomal proteins in early larval stages. This weak response is observed across diverse Wolbachia strains from both nematode and insect hosts suggesting a potential pan-Wolbachia transcriptional response during host development that diverged more than 700 million years ago.

Nearly Complete Genome Sequence of Brugia pahangi FR3.

Brugia pahangi is a zoonotic parasite that is closely related to human-infecting filarial nematodes. Here, we report the nearly complete genome of Brugia pahangi, including assemblies of four autosomes and an X chromosome, with only seven gaps. The Y chromosome is still not completely assembled.

Complete Genome Sequence of wBp, the Wolbachia Endosymbiont of Brugia pahangi FR3.

Lymphatic filariasis is a devastating disease caused by filarial nematode roundworms, which contain obligate Wolbachia endosymbionts. Here, we assembled the genome of wBp, the Wolbachia endosymbiont of the filarial nematode Brugia pahangi, from Illumina, Pacific Biosciences, and Oxford Nanopore data. The complete, circular genome is 1,072,967 bp.

Nearly Complete Genome Sequence of Brugia malayi Strain FR3.

Lymphatic filariasis affects ∼120 million people and can result in elephantiasis and hydrocele. Here, we report the nearly complete genome sequence of the best-studied causative agent of lymphatic filariasis, Brugia malayi The assembly contains four autosomes, an X chromosome, and only eight gaps but lacks a contiguous sequence for the known Y chromosome.

Sex chromosome evolution in parasitic nematodes of humans.

Sex determination mechanisms often differ even between related species yet the evolution of sex chromosomes remains poorly understood in all but a few model organisms. Some nematodes such as Caenorhabditis elegans have an XO sex determination system while others, such as the filarial parasite Brugia malayi, have an XY mechanism. We present a complete B. malayi genome assembly and define Nigon elements shared with C. elegans, which we then map to the genomes of other filarial species and more distantly related nematodes. We find a remarkable plasticity in sex chromosome evolution with several distinct cases of neo-X and neo-Y formation, X-added regions, and conversion of autosomes to sex chromosomes from which we propose a model of chromosome evolution across different nematode clades. The phylum Nematoda offers a new and innovative system for gaining a deeper understanding of sex chromosome evolution.

N6-methyladenosine regulates the stability of RNA:DNA hybrids in human cells.

R-loops are nucleic acid structures formed by an RNA:DNA hybrid and unpaired single-stranded DNA that represent a source of genomic instability in mammalian cells1-4. Here we show that N6-methyladenosine (m6A) modification, contributing to different aspects of messenger RNA metabolism5,6, is detectable on the majority of RNA:DNA hybrids in human pluripotent stem cells. We demonstrate that m6A-containing R-loops accumulate during G2/M and are depleted at G0/G1 phases of the cell cycle, and that the m6A reader promoting mRNA degradation, YTHDF2 (ref. 7), interacts with R-loop-enriched loci in dividing cells. Consequently, YTHDF2 knockout leads to increased R-loop levels, cell growth retardation and accumulation of γH2AX, a marker for DNA double-strand breaks, in mammalian cells. Our results suggest that m6A regulates accumulation of R-loops, implying a role for this modification in safeguarding genomic stability.

Author Correction: Large Enriched Fragment Targeted Sequencing (LEFT-SEQ) Applied to Capture of Wolbachia Genomes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Drug Repurposing of Bromodomain Inhibitors as Potential Novel Therapeutic Leads for Lymphatic Filariasis Guided by Multispecies Transcriptomics.

To better understand the transcriptomic interplay of organisms associated with lymphatic filariasis, we conducted multispecies transcriptome sequencing (RNA-Seq) on the filarial nematode Brugia malayi, its Wolbachia endosymbiont wBm, and its laboratory vector Aedes aegypti across the entire B. malayi life cycle. In wBm, transcription of the noncoding 6S RNA suggests that it may be a regulator of bacterial cell growth, as its transcript levels correlate with bacterial replication rates. For A. aegypti, the transcriptional response reflects the stress that B. malayi infection exerts on the mosquito with indicators of increased energy demand. In B. malayi, expression modules associated with adult female samples consistently contained an overrepresentation of genes involved in chromatin remodeling, such as the bromodomain-containing proteins. All bromodomain-containing proteins encoded by B. malayi were observed to be upregulated in the adult female, embryo, and microfilaria life stages, including 2 members of the bromodomain and extraterminal (BET) protein family. The BET inhibitor JQ1(+), originally developed as a cancer therapeutic, caused lethality of adult worms in vitro, suggesting that it may be a potential therapeutic that can be repurposed for treating lymphatic filariasis.IMPORTANCE The current treatment regimen for lymphatic filariasis is mostly microfilaricidal. In an effort to identify new drug candidates for lymphatic filariasis, we conducted a three-way transcriptomics/systems biology study of one of the causative agents of lymphatic filariasis, Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host Aedes aegypti at 16 distinct B. malayi life stages. B. malayi upregulates the expression of bromodomain-containing proteins in the adult female, embryo, and microfilaria stages. In vitro, we find that the existing cancer therapeutic JQ1(+), which is a bromodomain and extraterminal protein inhibitor, has adulticidal activity in B. malayi.